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Abstract 

The influence of particle-phonon coupling on pairing correlations in nuclei 
is studied by solving the Dyson equation including the anomalous (pairing) 
Green function. We develop the formalism for solving the equation with the 
minimum of approximations. The solution of the Dyson equation is compared 
with the diagonalization of particle-phonon coupled Hamiltonian in a small 
space. This comparison reveals that the effect of many-phonon states is incor- 
porated in the Dyson equation. We calculate the pairing gap of the neutron in 
120 Sn. We compare analytically the present method with a simpler treatment 
based on Bloch-Horowitz perturbation theory. 
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1 Introduction 



The study of the nature of pairing correlations is one of the most fundamental 
problems in nuclear physics. These correlations have been introduced in the study 
of nuclear structure in analogy with the case of metals U and since then played an 
important role in the understanding of many aspects of nuclear behaviour, because 
they affect static and dynamic properties of atomic nuclei, their structure as well as 
their reactions. One can quote, for example, the systematic staggering of the odd- 
even mass difference (section 2-1 in O), the reduction of the moment of inertia of 
collective rotational bands as compared to the rigid-body value (section 4-3 in [|J), 
the enhancement of the two-particle transfer reactions to specific final states |4j], the 
appearance of signals in both one-particle pick-up and stripping reactions for the 
same single-particle orbit near the Fermi level (section 5-3 in 0), the lifetimes of a 
decay in the heavy mass region || and so on. Influence of the pairing gap can be 
also found in the spectra of deformed nuclei. 1 Q Namely, the spacing of the low- lying 
non-collective excited states (band heads) in odd nuclei is reduced in average due to 
the pairing gap ||. Two-quasiparticle states are of relatively high energy in heavy 
even-even nuclei, because breaking the pair costs energy. The feature analogous to 
the latter can be seen also in odd nuclei (sections 5-2 and 5-3 in 0). The interplay 
between the pairing and the other degrees of freedom also plays an important role. 
For example, pairing correlations are suppressed by rotation ]7|. § . 

Since pairing correlations play such an important role, it is crucial to ask what 
is their origin. 

Pairing correlations in finite nuclei have mostly been studied in mean field the- 
ory through BCS and Hartree-Fock-Bogoliubov (HFB) calculations, making use of 
phenomenological interactions, the monopole pairing interaction having probably 
been the most often used (e. g. 0). The list of the interactions may include the 
quadrupole pairing interaction (used for improving the monopole pairing interac- 
tion, see e. g. []T(J for rotation and |Tl[ for the pairing vibration.), the surface-delta 



interaction e. g. [12|, (13] and [14j], density-dependent zero-range pairing forces [15 
the Gogny force, |Tf| and (T7J, and so on. 

Recently, several studies have been devoted to pairing in nuclei near the drip line, 
for which it is crucial that the pairing force treats correctly the interplay between 
the discrete and continuum unperturbed single-particle levels correctly: that is, 
summation of the unbound components of many-body wavefunction caused by the 
interplay must be canceled out |Tj|,[jrjJ. 

The connections of pairing correlations in finite nuclei with the bare interactions, 
that is, on interactions which reproduce the properties of scattering of free nucleons, 
have been very scarce in finite nuclei (for a few exceptions, cf. 



20 , 21 and 12 



On the other hand, this represents an important line of research in nuclear matter 
and references therein 



27 



and in its applications to nuclear astrophysics, 
where for example pairing plays a very important role in the cooling of the neutron 
stars by neutrino emission (e. g. [ffifl ). Also much studied have been the effects of 
medium polarization on pairing correlations. In terms of field theory this effect is 
given by bubble diagrams of the particle-hole excitations p9|-p4|. 

1 The influence can be recognized in average for low-lying levels. Evidence of the correlations 
in each level is not always easy to find due to state dependence of the pairing gap. 
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In the limit of infinite number of the bubbles, one can think of an interaction be- 
tween nucleons mediated by the exchange of phonons (phonon-induced interaction) . 
In this idea the coupling, or vertex, between particle 2 and phonon is crucially 
important, and the influences of this coupling on many properties of phonons as 
well as of particles have been investigated intensively in finite nuclei, although not 
its effect on pairing correlations. In particular, the coupling with low-lying surface 
vibrations renormalizes the single-particle motion in an essential way, shifting the 
energy of the levels and changing the level density and the effective mass around the 
Fermi surface, and producing a significant breaking of the single-particle strength 
for the levels lying far from the Fermi surface ||30|| , p5 1 . 



A formulation which can properly investigate the effect of the particle-phonon 
coupling on pairing correlations in nuclei was already proposed in 1960 's by Belyaev 



et al. |36], [37], [3g] and Migdal et al. p9] , f40f[ on the basis of field theory. The key 
point of the method is to introduce anomalous (pairing) Green functions in addi- 
tion to normal Green functions. In condensed matter physics this formulation was 



introduced by Gor'kov |^T| and Nambu [12], and the original dynamical equation 
was rearranged with approximations characteristic in the metal and is known as 
Eliashberg equation or Gor'kov equation, which was shown to be powerful in 
the calculation of the tunneling probability. A relatively recent review is given by 
[H. The application of the particle-phonon coupling to pairing problems in finite 
nuclei has been, however, surprisingly rare. Recently Barranco et al. solved BCS 
equation with the phonon-induced interaction obtained according to Bloch-Horowitz 
perturbation ES| . They have also shown the crucial importance of this coupling for 
the stability of the halo nucleus 11 Li [46 1. Avdeenkov and Kamerdzhiev |47j solve 



the Dyson equation in the linear approximation. ( See also [Q . ) Keeping in mind 
the historical background described above, we investigate in this paper the effect of 
the phonon-induced interaction on nuclear pairing correlations. 

We solve the Dyson equation consisting of the one-body as well as anomalous 
Green functions for particle-phonon coupled systems and calculate the pairing gaps. 
The phonon is treated as an independent degree of freedom. In nuclei, this is a tech- 
nique to treat complicated diagrams of interaction nonperturbatively. Introduction 
of the anomalous Green function can be understood also as this kind of technique. 
It is emphasized that in the field-theoretical approach it is possible to include effects 
beyond the mean field as well as the mean pair-field approximation, which mani- 
fest themself in the fragmentation of the single-particle strength. After making a 
choice of proper self-energy diagrams, we solve Dyson equation with the minimum 
of approximations, being careful not to suppress the effects beyond BCS. Thus the 
most significant point of this study is that we investigate the pairing correlation 
without assuming a priori that the single-particle or quasiparticle picture is good. 
In addition a class of diagrams is included nonperturbatively in the calculation. We 
shall pay attention also to relation between this approach and diagonalization of 
particle-phonon coupled Hamiltonian in the many-body plus phonon space. 

Finally let us mention two more groups concerning the Dyson equation in nuclear- 
structure. One of them is Muther and Dickhoff et al. They developed an approx- 



2 This terminology is used in different ways depending on the context. One is a single-particle 
sitting on a level above the Fermi level, and another is simply single-particle or nucleoli. In this 
paper we use the terminology mainly in the latter sense. 
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imate way to solve Dyson equation ||9j ( without the anomalous Green function ) 
and applied their method to doubly-magic nuclei |50fl-|]54j| with a nucleon-nucleon in- 



teraction, emphasizing fragmentation of the single-particle strength. An application 



was made also to nuclear matter |55 |. Another group which we note is Waroquier 
et al. They proposed another way to approximate a solution of Dyson equation 
and applied it to some nuclei [0, |57|, ^8|. They performed calculation also for an 



open-shell nucleus |59[ using a formulation which includes a concept of the pair- 
ing gap, however, the gap was not discussed in their paper. This group uses both 
nucleon-nucleon interaction and particle-phonon coupling. A comparison was made 
of the two methods of the groups in fl60fl . 

This paper is organized as follows: In section 2 we discuss the formulation and 
describe our way to solve the Dyson equation. We also consider the relation between 
diagrams with and without the anomalous Green functions. In the next section 
detailed comparisons are made between the solution of Dyson equation and the 
diagonalization of the particle-phonon coupled Hamiltonian. Section 4 shows results 
of the calculation of the neutrons in 120 Sn using Dyson equation. We discuss the 
spectral functions, pairing gaps, the single-particle energy shifts and the pairing 
energy. BCS approximation is discussed briefly in section 5. The last section is 
devoted to a summary. The equations used in this paper are derived thoroughly 
in appendices. Appendix A is for the derivation of the diagram rule. The self- 
energy used in the calculation is derived in Appendix B. The next appendix is 
about the general form of the particle Green functions, and some useful relations 
are derived. Appendix D treats equations for calculating the poles and the residues. 
The equations of the total energy and the pairing gaps are derived in Appendices E 
and F, respectively. 

2 Formulation 

The core of this study is the Dyson equation 

g; 1 M = gJ- 1 m-s,m, (i) 

where each term is a 2x2 matrix, i. e., 

Gu{uj) = 1 , (2) 

G >) = 21 22 h ( 3 ) 

E>) = ( W W j ' (4) 

with uj being an energy variable. G>(co>) is a perturbed nucleon Green function 
obtained by solving the equation, fi denotes the spherical good quantum numbers 
(nlj)^ of a single-particle orbit throughout this paper except for Appendix A which 
treats the most general case. We assume that the system considered has spherical 



4 




(I (I 



Figure 1. The self-energy, uj and uj' denote the energy variable. A stands for 
the multipolarity of phonon, and /j, and fj,' are the single-particle indices. 



symmetry. The off-diagonal elements of G^(uj) are anomalous Green functions. 
G° {yj) is an unperturbed nucleon Green function 



1 



(5) 
(6) 
(7) 



where e° denotes an unperturbed single-particle energy, and Ef is Fermi level. r/^ is 
a parameter as 

where 77 is a real positive parameter. In this study we consider a self-energy X^(cu) 
given by the diagram in Fig. 1. 

The wavy line in Fig. 1 denotes an unperturbed phonon Green function 



h 



+ 



uj — uj' — hQ\„ + ir)D —oj + uj' - hVt\ n + ir\ D 



(9) 



where hfl\ n is the phonon energy with the multipolarity A and is uniquely specified 
by another label n, if many modes have the same A. t]d stands for a real positive 
parameter. The thick line indicates the perturbed nucleon Green function |Gy(u/). 
The small open circle is an unperturbed vertex of the phonon and the nucleon: 



(") 



(10) 



B\ n is the inertia parameter of the phonon. R and U (r) are a nuclear radius and 
a nuclear potential, respectively. From the diagram we obtain the equation of the 
self-energy 
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h 



l 



X 



p 

(/ill^o^nil/i') 



2fi Ari E Ari 2^ + 1 



ft 



(11) 
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where t% = ( J -1 ) ■ The derivation of the diagram rule is given in Appendix A, 
and the self-energy in a spherical nucleus is derived in Appendix B. 
For solving Eq. ([I]) we put 
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aa 
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Assuming the time-reversal invariance of the ground state, we can put 



^22 




-iu)rj 



(12) 
(13) 

(14) 
(15) 



where uj G a + denotes a pole of the nucleon Green function ( Reuj G a + > ). a is a label 
to distinguish poles associated with /i. Rj} a (±uj G a + ) are residues of Gj}(uj)/h at the 
poles ±uj g \. Eqs. (|T2|)-(|T5|) are explained in detail in Appendix C. Inserting the 
above equations to Eq. flTT|) we obtain 

11 <^i^w> 2 1 



ill, 



EE 

fj,'a An 



2Q\ n B\ n 



dr 



2J, 



1 



-co: 



G+J 



EE 

u'a An 



+ 



(16) 



(n\\Ro^Y x \\^) 



-1 



+ 



2j u + 1 



An 



UJ 



UJ. 



fi'a 
G+ 



fia 



(17) 

ir]D ' 
's are real, then 



S 22 ^) and TP^{uj) can be obtained in the same way. ( If i?* 2 (±u; ( 
Si 2 (a;) = Ej 2 (-a;). ) 

Keeping the idea that stationary states are treated, we calculate u G a + and R^ a (±uJ G a + ) 
etc. in the following way. First ±Re uj G a + are determined by searching for the points, 
numerically on the real uj axis, satisfying 

detG~\±ReuJ G a + ) = , (18) 



where 



with 



hc;\uj) 



uj-e^-Reh^iuj) 
~^f(uj) 



e° 



-hz^iuj} 



ReKLf(uj) 



(19) 



e ^ — ep. The above equation is the same as UG^i^uj) except for the non- 
hermitian component of the diagonal elements eliminated. We identify the residue 



of G l hu)/n with that of G]^(uj)/h: 



R aai.^ Ujl G+ 



uj + e® -hReZf(uj) 



-fh 2 det G„ ! 



(UJ) 



(20) 
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where the derivative is taken on the real axis. For the pairing part the residue is 
calculated 

R™W- 



(21) 



The imaginary part of the pole may be approximated by 



Im^al ~ Im __A_J$(a$.) , ( 22) 



or one can set Im 0Jq\ equal to a small constant. In section 4 the latter method is 



used. Eqs. (|2Xf )-(f22|) are an extension of the argument in the sections 3.4 and 4.3.6 
in |35|], and are discussed in Appendix D. 

Eq. (Q) can be solved by iteration using the BCS Green function ( section 7-2 in 
[[£J ) as an initial guess of a solution: 

IcfHou) = ( $ + $ )e ir >», (23) 

(2 2 \ 
^ + ^ ) e"^ , (24) 

\ G f\u) = - U /r + U ;\. , (25) 

k * UJ + £ M — IT] iO — + IT] 



where £ M = y (e° — ef) 2 + A 2 , and v M being the occupation factors of BCS theory. 
A is a parameter. The BCS Green function is inserted to Eqs. (|I~6|) and (|TTD, i.e. 



S^Re^ + , (26) 

— r] — > Imug + , (27) 

vl^Rf(-^ G+ ), (28) 

uj-<(wg+), (29) 

W - . (30) 

Note that Eqs. ([T^) and flT7] ) do not have the summation J2 a a ^ this stage. Then 
the new G^uj) can be calculated by the matrix inversion of the right-hand side of 
Eq. ([IT). Finding the real points satisfying Eq. (|TSJ) and calculating Eqs. ( PU[ ) and 
(pT|), one can get new poles u}q" + and .R^io;^), (2, j = 1,2). These are substituted 
to Eqs. (|i~6l) and (|T7|) again, and the procedure is repeated. The converged G^u) 
and S At (cu) give a solution of Eq. (Jl|). is determined so as to obtain a right 
expectation value of the nucleon number (N) for the nucleus under consideration. 

It is worthy to discuss what would happen if the energy denominators in Eqs. (|16| ) 
and fllTD become very small around the poles of the particle Green function. This 
can happen when the unperturbed energy levels are very dense. Let us ignore pairing 
correlations for simplicity, and consider the pole -w^. Then we have 



1 - 



(31) 
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Thus if the energy denominator of (no) is very small around the pole, then so is 



^m(~ 0J G+)- On the other hand the sum rule of the single-particle strength must be 
satisfied: 

E(^K" + ) + «)) = i. (32) 

a 

Therefore if the small-energy denominators are encountered so often, then the single- 
particle strength becomes distributed among many poles, leading to a strong frag- 
mentation [fjOfl . This mechanism means a technical advantage of Dyson method 
because the small-energy denominators do not cause severe problems as can happen 
in Rayleigh-Schrodinger perturbation theory. 

If the a;-dependence of the self-energy £ M (cu) is ignored, then detG^ (a;), see 
Eq. (|I9|) , is a parabola, and there are only two roots in Eq. (|18|) — that is BCS. With 
the (^-dependence of £ M (w) switched on, det (u) behaves in the way schematically 
shown in Fig. 2a. 

It is seen that the singularities of S^(o;) create many poles of G^{uj). In Fig. 2b 
we illustrate a magnification in a neighborhood of a singularity. When the averaging 
parameter t] and t]e> are finite, the determinant is a continuous function on the real 
to axis. It is noted that the zero-point A in Fig. 2b is a spurious pole, because 
the curve becomes discontinuous at the point A with i],i]e> —>■ 0, and thus the pole 
disappears. In our experience the residue R^^ujq^) for the point A, calculated 
using Eq. (p0[), is usually negative. 

We note that effect of the off-diagonal Green functions G(nij)(n'^nij) are assumed 



small and neglected in this paper. ( Cf. |59| and p3| 



Now let us consider a field theoretical question related to the anomalous Green 
function. Is it possible to express diagrams including the anomalous Green functions 
in terms of the normal Green functions? Let us assume an equation 





(33) 

The left diagram is a two-body Green function carrying the pairing correlations 
expressed without the anomalous Green functions (see section 34, chapter 7 in []61~| ). 
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The ladder diagram of the phonons may be included in the complicated diagram. 
The right-hand side indicates the anomalous Green functions in the notation of 
Migdal. ( Cf section 1.3.4 in [39|. ) In this discussion, a line represents not a 2 
by 2 matrix but either a normal or a anomalous nucleon Green function. One of 
components of the correlation energy associated with the solution of Dyson equation 
can be written, up to a factor, as 




> 








1 




- h- 




) 



(34) 



where the thick lines are perturbed Green functions, and the thin lines are unper- 
turbed ones. Therefore it is seen that the diagrams of the following type are included 
in the proper self-energy in the Dyson equation: 




(35) 



See also chap. 2 in 



3 Comparison with a diagonalization of the particle- 
phonon coupled Hamiltonian 

In order to understand this field-theoretical calculation from a viewpoint of a method 
more familiar to the nuclear-structure physics, we made a comparison with a diag- 
onalization of the Hamiltonian 

H = ifp ar + H coup + ifpho , (36) 

where if par and if p ho are unperturbed particle and phonon Hamiltonians, respec- 
tively. The particle-phonon coupling H coup was taken from H' in section 6-5a in 
[H]. We diagonalized if in a space of two particle ® 0-2 phonons coupled to the 
angular momentum J = 0. For the particle space we took { g 9 / 2 , d 5 / 2 , J15/2 } above 
the N = 126 shell gap calculated with a Woods-Saxon potential ]4"5| . For phonons 



we used the lowest 2 + and 3~ solutions of a QRPA calculation performed using the 
multipole-multipole force with strengths adjusted so as to reproduce the observed 
transition probabilities in 208 Pb. The coupling strength can also be obtained from 
the QRPA calculation (cf. sections 6-2c and 6-5a in ||). 

The input energies and vertex are shown in Table I and II, respectively. 
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Table I. The input energies used in the calculation, is the unperturbed 
single-particle energy, and h£l\ is the phonon energy. 



fi 


[MeV] 




TiVLx 
[MeV] 


g9/2 


-4.314 


2+ 


4.10 


jl5/2 


-2.490 


3" 


2.10 


d 5 /2 


-1.992 







Table II. The matrix elements of the vertex. The vertex matrix is real and 
Hermite. 



n 




A" 


[MeV] 




g9/2 


g9/2 


2+ 


1.36 


d 5 /2 


d 5 /2 


2+ 


0.87 


jl5/2 


jl5/2 


2+ 


2.24 


g9/2 


d 5 /2 


2+ 


1.26 


jl5/2 


g9/2 


3" 


4.02 



We compare two kinds of quantities. One is the occupation probability of the 
orbits 

v\ (Dyson) = ]T(c 

= E4i(-^ a + ), (37) 

a 

<(Diag) = EiVxl^ + EK^^T. (38) 



where a M 2 X and a M y are the amplitudes of the components of the ground state |gs) 
in the diagonalization method: 

\gs) = J2^x\X) + J2^v\Y) , (39) 

X Y 

\X) = [|^ 2 >® |0, lor 2 phonon)] j7r=0+ , (40) 
\Y) = [\fi)® \n'±n) ® |0,lor2phonon)]j^ =0 + . (41) 

Another quantity to compare is the correlation energy 

E cor = E - E° unp , (42) 
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Table III. The occupation probabilities of the single-particle orbits v 2 and the 
correlation energies E COT of the two methods, j] = r/D = 1 keV and £f = 
—4.7908 MeV were used in the Dyson calculation. 







v 2 




E 

-^cor 




g9/2 


jl5/2 


d 5 /2 


[MeV] 


Dyson 


0.178 


0.013 


0.001 


-1.62 


Diagonalization 


0.179 


0.013 


0.001 


-1.53 



Table IV. Another set of the input energies. 



n 


[MeV] 




[MeV] 


g9/2 


-3.00 


2+ 


0.50 


jl5/2 


-2.50 


3" 


0.50 


d 5 / 2 


-2.00 







where E® np is the unperturbed ground-state energy, which is equal to 2e° 9//2 in the 
present model. The total energy E is given in the Dyson method by 

E = Y,Rll(-^ + )(^ + l)(-^c + )+e F (N) . (43) 

fia 

This equation is derived in Appendix E. The counterpart in the diagonalization 
method is 

E = the lowest eigenvalue of H — (H pho ) . (44) 

The result is shown in Table III. 

It is seen that the two calculations give a nearly-identical result. In order to 
see difference between the two methods, we used another set of the unperturbed 
energies artificially changed ( see Table IV ) keeping the matrix elements of the 
particle-phonon coupling unchanged. 

This choice was made in accordance with the discussion on the fragmentation in 
section 2. The physical difference between the two models is clear in the distribution 
of the single-particle strength in Fig. 3. 

Given that there are only three particle levels, the latter model (Table IV) is an 
extremely fragmentation-enhanced model. The correlation energies of the second 
model is plotted in Fig. 4 (p. 13). 

We emphasize that Dyson equation is solved nonperturbatively including many 
poles, so that the many-phonon diagrams are included in the solution. Figure 4 
clearly shows this many-phonon effect. It is also noted that the present solution 
of Dyson equation does not have the vertex correction, while the diagonalization 
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Figure 3. The single-particle strength distribution of the three orbits. u)q" 
means ±Wg a + . The left three panels are results of the first model (Table I), 
and the right three are of the second model (Table IV). Note the difference in 
the vertical scale. 



method includes the effect within the order of 2-phonon. Thus the truncation scheme 
in the actual calculation is in fact different between the two methods. Given this 
difference, the comparisons show that the two methods are consistent. 

In our solution, the number of the poles of the particle Green functions can 
vary from one step to another in the iteration process, and due to this reason small 
fluctuations are unavoidable. In most calculations presented in this paper, including 
those of the next section, the errors in the pairing gaps and perturbed single-particle 
energies of the quasiparticle poles in the valence shell are negligible. 

4 Pairing gaps in 120 Sn 

In the wake of the above calculation, we have performed a more realistic study, 
calculating the neutron pairing gap of 120 Sn using the Dyson equation. (A part of 
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Figure 4. The correlation energies of the Dyson equation ( E cor = —3.55 =t 
0.01 MeV, the mean value is shown by the line. For the error, see text. ) and 
the diagonalization method (the crosses) obtained for the input in Table IV. 
The horizontal axis is the maximum phonon number in the basis states used 
in the diagonalization, e. g., phonon number = 1 means that the basis state 
is [ \nfj,') (8> |0or lphonon) ] j 7r= o+. E COT of the Dyson method corresponds to a 
large phonon number (see the text). £p = —4.280 MeV was used, rj and ijd 
are unchanged. 



this calculation has been published in a compact way in [64]. ) In this calculation the 
single-particle basis covers all of the bound energy region starting from the orbital 
lsi/2- We used the unperturbed single-particle spectra 

el=— (e™-e F )+e F , (45) 
p m k p 

with mt = 0.7m and e^ s being Woods-Saxon spectra, for avoiding double counting 
of the particle-phonon coupling effect. ( Cf. section 4.6.3 in [)55| . ) 

The computation time needed for our calculations depends strongly on the num- 
ber of phonon modes An included. The full QRPA response for the multipolarities 
A 71 " = 2 + ,3~,4 + ,5~ in the energy interval 0-20 MeV used in ref. |45[ consists of 



about two hundreds phonon modes of energies TiVt\ n and zero-point amplitudes f3\ n 
for each multipolarity. We include the four lowest phonons, one for each multipo- 
larity, which give the largest contributions to the induced phonon interaction. We 
account for the effects of the other roots including only a few effective phonons of 
energy M2^, distributed in the interval 0-20 MeV, choosing their effective strength 
so that when they are used in the Bloch-Horowitz calculation of ref. |^5[ they re- 
produce the state-dependent gap obtained there. This is obtained, considering that 
the sum of the (asymmetrized) matrix elements of the induced interaction between 
two pairs (j u )j=o, (ji/)j=o due to the phonons lying in an energy interval [f2 a ,0&], 
calculated according to the Bloch-Horowitz formalism, is given by 

IWjgf Hi 2 £ 



V — — (46) 

where E cor and e u are the correlation energy of the ground state and the absolute 
value of the single-particle energy with respect to the Fermi level, respectively. The 
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Table V. The energies of the phonon modes frfi^ and coupling strength 
(3fj^/2\ + l. Tables a), b), c) and d) are for A 71 " = 2 + , 3~, 4+ and 5~, 
respectively. The lowest-energy modes (n = 1) were taken from a QRPA 
calculation directly. The coupling strengths as well as the energies of 
the other modes were determined by the procedure shown in Eqs. ( f46[) 
and (47). 

a) b) 



n 



A 71 " = 2+ A 71 " = 3" 

[MeV] [MeV] 



1 1.173 0.0554 1 2.423 0.0591 

2 5.2 0.0134 2 5.57 0.0317 

3 12.5 0.0447 3 10.0 0.0238 
4 21.0 0.0291 



c) d) 



n 


A 71 " 
[MeV] 


= 4+ 

0f /V2A + 1 


n 


A 71 " 

hO cS 

nlL \n 

[MeV] 


= 5" 

(3f/V2X + 1 


1 


2.470 


0.0248 


1 


2.402 


0.0250 


2 


8.0 


0.0300 


2 


8.0 


0.0365 


3 


12.0 


0.0300 


3 


13.0 


0.0166 


4 


15.0 


0.0270 


4 


21.0 


0.0232 



effective strength of the phonon representing this interval is then chosen so as to 
satisfy the equations 

(f)eS\2 n b n2 



E COI - {e v + e v > + httfj n x jtn a E ™* ~ ( e » + + h£l\ nj 

hnf n = hQ b , (47b) 

for the pairs (j u )j=oi Cv)j=o giving the largest contribution to the pairing gap for 
the multipolarity A. The energies and zero-point amplitudes of the effective phonons 
are listed in Table V, divided by v / 2A~+T. 
Then the coupling strength 



is used for Eq. (0). The BCS+Bloch-Horowitz calculation performed with this 
restricted ensemble of phonons reproduces the state-dependent pairing gaps of ref. 



45] within a few per cent. 
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Below, results obtained with the fixed imaginary parameter f] = rjQ = -Imw^ = 
1 keV are shown. In the present calculation the maximum number of poles with re- 
spect to n is ~200. The sum rule of the single-particle strength, Eq. (0), is satisfied 
by more than 90 % in average in the orbits in the valence shell. The computations 
have been performed on a parallel computer. In fact the Green function method is 
well suited for parallel computation, because the pole search of each orbit can be 
done separately. 

A few self-checks are possible on the accuracy of the solution of Dyson equation. 
The self-energy is a functional of the particle Green function, and if this Green 
function, represented by Eqs. (|l2"|) -( |l5l) , is identical to {G° (u) — X M (co>) )~\ then 
the Green function is a solution. Such a direct check is shown in Fig. 5 for the orbits 
in the valence shell around the Fermi level. 

It is seen that the accuracy is good. Another check is Eq. ( |136| ) in Appendix C. 
This equation is satisfied with good accuracy. 

We show the calculated pairing gaps in Fig. 6 (p. 17). 

The value of the gap of 2ds/2 is mostly determined by the off-diagonal element 
of Eq. (0) with n' = lhn/2 and A 71 " = 3~ Therefore \R{l n/2 (ojG+ /2a *)\ = 0.35 of 
the quasiparticle pole 3 f\ a plays an important role not only for A lhll / 2 but also for 

^2d5/2- On the other hand |^2d5/2( w G+^ 2a °)l * s 0.04, thus the influence of 2d 5 / 2 on 
A lhll / 2 is small. 

The average of the pairing gaps in the valence shell is 0.54 MeV in the present 
calculation. This may be compared with the value of 1.39 MeV deduced from the 
observed odd-even mass difference as well as the average gap of 0.82 MeV obtained 
in ref . . 

The poles ±^4- °f the Green function G 11 are distributed symmetrically with 
respect to the Fermi surface (uj = 0). The spectral function of the orbits in the 
valence shell are illustrated in Fig. 7 (p. 18). The single-particle picture is rather 
good for most orbits, except for the d$/2 whose sinlge-particle strength is split into 
several peaks. Most orbits have a pronounced hole- or particle-character, except 
for for the /in/2 which displays two large peaks on either side of the Fermi energy. 
It is to be noted that for this orbit the single-particle strength for the two quasi- 
particle poles is only 0.73. The third not negligible peak close to uj =2 MeV may 
be associated with the quasiparticle plus a collective state. 

It is seen from Eq. (|1 12|) in Appendix C that the peaks in uj < part of the 
spectral function are associated with a pick-up reaction and those in the uj > part 
with a stripping reaction. Thus the spectral functions of Fig. 7, can be rearranged 
as shown into Fig. 8 (p. 19) for the case of h u / 2 . 

If each excited state can be resolved, the definition of E x is Eq — Ei, where 
Eq is the ground-state total energy of the neutrons in 120 Sn relative to sf(N), and 
Ei is defined in the same way as Eq but for various states of 119 Sn and 121 Sn. If 
strong signals are observed in both pick-up and stripping reactions for the same 
orbit near the Fermi level with the same E x , this would be a clear indication of 
pairing correlations. 

The particle-phonon coupling influences of course also the single-particle spec- 



3 A pole, or a pair of poles if the pairing gap is not zero, which carries the largest, or major, 
single-particle strength R]} a (uj^ a + ) + R]} a {-uj^ a + ). 
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Figure 5. Absolute value of the real part of Gj}(uj)/h in log scale. The left 
column shows 
those given by 
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, and the right column illustrates 
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eg [MeV] e^ 1 [MeV] 

Figure 6. a) Absolute values of the pairing gaps of the neutron of 120 Sn in 
the calculation of the Dyson method shown as a function of e°. The gap 
is obtained according to Eq. ( |158| ) in Appendix F. The arrow indicates the 
location of the Fermi level, b) |A^| as a function of the perturbed single- 
particle energy e\ (see Eq. (|160| ) in Appendix F) in the vicinity of the Fermi 
level. 



trum. We show the perturbed as well as unperturbed spectra in Fig. 9 (p. 19). 

The increase in the level density in the vicinity of the Fermi level is a known 
effect of the coupling (section 4 in pp|). 

We calculated also the pairing energy (see Appendix E): 

(48) 

and obtained -E n air = —3.9 MeV. 



5 BCS approximation 

This approximation was introduced in section 2 as the initial guess for solving Dyson 
equation. From the viewpoint of Green function method, BCS approximation is a 
normalized-quasiparticle approximation. Thus, generally speaking, if the single- 
particle strength is well concentrated on the quasiparticle pole, then BCS may be 
a good approximation. As was discussed in section 2, if the self-energy has no 
a;-dependence, then the single-particle strength is distributed to only two poles at 
most. 

Next we compare analytically the present Dyson calculation and the BCS + 
Bloch- Horowitz calculation in |fI5| |. One of the differences is, as was mentioned, that 



Rj^, (^G+) (do denotes the quasiparticle pole) differs from u^v^ in BCS due to the 
fragmentation of the single-particle strength. The second difference is the deviation 
of the function Z^Uq 1 ?), Eq. ( |154|) in Appendix F, from 1. A third difference is the 
energy denominator in the gap equation. If one puts for simplicity TlT}^{uj^) ~ A M , 
Ruao( u G+) — u n v n an d ^67+ — *he pairing gap in Dyson method can be written 
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Figure 7. The spectral function S^lo) = Im G 1 ^{uj)/Tl /it for the orbits in 
the valence shell. The poles and residues of the solution were inserted into 
Eq. (12) with Imw^ = —100 keV. u = corresponds to the Fermi level. 
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On the other hand, the gap of BCS + Bloch- Horowitz calculation leads to 

h l 



(49) 
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Figure 8. The spectral functions of the neutron as functions of E x 

a) hole type b) particle type. 
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Figure 9. The unperturbed (e°) and perturbed (e*) single-particle spectra 
associated with the quasiparticle poles, is defined in Eq. (|160|) in Appendix 
F. 



^co r -(|^| + |^|+^An) ' ' ' 

where _E cor is the correlation energy of the perturbed ground state. Thus, the es- 
sential difference, in the no correlation limit, is the sign of in the second energy 
denominator in Eq. (|4"P|). 

Eq. (fl8|) may be compared to the pairing energy in the BCS approximation: 

£ pair (BCS) = }^ — ~ — u^VV^V ' ( 51 ) 

VftH' = ((^ / ^0(/ i/m M')l^ / |(A im M)(/ im M))antisymmctri z ed • (52) 

We obtained £ pair (BCS) = -4.6 MeV in the BCS + Bloch- Horowitz calculation. 



6 Summary 

In this paper we have investigated the effect of particle-phonon coupling on nuclear 
pairing correlations in detail. The dynamical equation treated is Dyson equation 
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including the anomalous Green function. We showed the formulation from the view- 
point to clarify how to solve the equation in such a way that the effect beyond BCS is 
not suppressed. The comparison was made between the solutions of Dyson equation 
and diagonalization of the particle-phonon coupled Hamiltonian, and we clarified 
both difference and similarity. The difference is in the truncation scheme, and non- 
perturbative effect can be taken into account for a class of diagrams in Dyson method 
by solving the equation iteratively. The similarity of the solutions is clear when the 
model does not have strong the fragmentation effect of the single-particle strength. 
We have solved Dyson equation for the neutrons in 120 Sn and calculated the pairing 
gaps. The average pairing gap near the Fermi level arising from the phono n-induced 
interaction is around 40 % of the observed gap from the odd-even mass difference. 
According to the spectral functions the quasiparticle picture is not perfect but ac- 
ceptable for many of the single-particle orbits near the Fermi level except for 2ds/2. 
This level is more than 2 MeV below the Fermi level, however, the orbit has the ap- 
preciable magnitude of the pairing gap. This gap arises mainly from the 3~ phonon 
coupling, of which the vertex matrix has only off-diagonal elements. For this reason 
the pairing gap as a function of the unperturbed single-particle energy does not have 
a completely-smooth shape peaked at the Fermi level. As for the overall shape, the 
pairing gap arising from the phonon-induced interaction is of the surface type. 

We have discussed also the BCS + Bloch-Horowitz calculation. The mechanism 
to avoid divergence possibly arising from small energy denominator is different be- 
tween Dyson calculation and BCS + Bloch-Horowitz method. The latter method 
has always negative, i. e. non-zero, energy denominators in the pairing gap. 

Several further steps are needed, in order to obtain a complete quantitative 
assessment of the contribution of the particle-phonon coupling to the nuclear pairing 
correlations. First, the vertex correction of the self-energy should be investigated. In 
fact, our present calculation is partially nonperturbative and partially perturbative. 
Needless to say, it is a nonperturbative effect that the pairing gaps were obtained 
using the single-particle basis. On the other hand, the vertex was treated in the 
lowest perturbation. 

It will be necessary also to combine the phonon-induced interaction with a bare 
NN interaction into the self-energy for understanding the nuclear pairing correlations 
fully microscopically, and to look for an adequate Hartree-Fock basis consistent 
with the effective mass arising from the particle-phonon coupling. One may be also 
concerned with breaking of the particle-number conservation. Since we calculate the 
pairing gap in a finite-body system, the particle number is not conserved. 



Appendix A — diagram rule — 

The original Nambu-Gor'kov formulation assumes that the spin up or down is a 
good quantum number [43 . In this appendix we show in detail a derivation of the 
diagram rule without this condition. We derive the rule by making an expansion of 
a perturbed nucleon Green function using products of unperturbed nucleon Green 
functions [p^j . It is noted that the nonperturbative effect of Dyson equation is 
essential for getting the anomalous Green function if unpaired states are used as the 
basis. 



20 



We put 



= ( J 1 ■ ( 53 ) 



In this appendix v is a single-particle index of a complete basis, and V is a time- 
reversed state of v. Note that {z/j} i= ... = {Ui} i= ..., and we do not impose any explicit 
good quantum number on the basis as long as a general diagram rule is discussed. 
Let us consider the particle-phonon coupling Hamiltonian ( see sections 6.2b and 
6.3a in § ): 

Hint = E «I«Im E Wt>>^73** > ( 54 ) 



where 



1 dU(r) 

7 im = ~Ro -^Y lm (9,p) , (55) 

a im = (ai)o (cf m + c^) , (56) 



( ai ) = <H^ ra |o> = y2^' (57) 



drr R — — Rq— - — . (5£ 



-R is a mean nuclear radius, and U(r) is a nuclear potential. c\ m is a creation 
operator of a phonon specified by the multipolarity / and its ^-component m. We 
define = (— ) m q_ m . (cty) is a zero-point amplitude of the phonon, which is 
related to the inertia parameter B\ and the phonon energy TiVLi. \lm) and |0) are a 
one-phonon and zero-phonon state, respectively. po( r ) denotes a nuclear density in 
the equilibrium. 
With a relation 

</i|F/» = (77|F/J7J>, (59) 

it follows that 

E^l^lk)^^^ = ^\FLW)kic u - CjI 4) (60) 

= Y,(^\ F LW)clc u + const . (61) 

That is, since the particle-phonon coupling considered here does not depend on the 
spin, Nambu-Gor'kov formalism can be introduced without much extension. 
We introduce a time-dependent coupling in the interaction picture: 

HUt) = E^^m(0E^l^lk)|*i(^3^) , (62) 

where the time-dependent operators are defined as 

ai m (t) = e< tln a lm e- ltl ^ , ( 63 ) 
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Hq 11 and Hq are unperturbed Hamiltonians of the phonon and single-particle relative 
to the Fermi level, respectively. The basic formula from which we start is 

iG(t) = 



iG^{t) iG%{t) 



m=0 



1 



m\ J-oo 



dti 



dt r 



X($ |T [#int(*l)#int(*2) " ' ' H iat (t m )*„(*) *+] | $ )connected • (65) 

T[- • •] is the time-ordered product, and the meaning of "connected" shall be clearest 
in the diagram later. |$ ) denotes an unperturbed ground state of the system of the 
single-particles and phonons. It is assumed c M |$o) = Qml^o) = 0. The m = term 
is ($o|T[\l/y(t)\l/J t ] |$o) = iGo(t). Since the unperturbed Hamiltonian does not have 
the coupling, we have 

l$o) = |^)K), (66) 

where and |$g h ) are the ground states of the single-particles and phonons, 
respectively. The derivation of the formula (^) is not affected by the inclusion of 
the anomalous Green function. (See, e. g., |)2] for the derivation.) 

By inserting Eq. (|62|) to Eq. (^) and using Wick theorem, it follows that 

iGll(t) iGg(f) 
iGllit) iG*(f) 

CO 

E 



m=0 



n 




dti 



dt m „ 
— oo ^ 



X E K h E^ll^imJ^l) 



E K lm E (PnWl 



t 



If-nTYl 



rn i ' 1 m 



x($o n |T [dj imi (ti)Qf fatn9 (t a ) ■ • ■oii mmm (t m )] m h ) 



x($S|T 



*I 1 (*l)T3*, 1 (tl)^ 2 (t 2 )T3^ 2 (t2) ■ • ■*l n (t m )T 3 *» m (t m )V v (t)*l \^) 



E 



l 

m! 



dt r . 



m=0 

x E E^il-^lmiK) ••■ E «<m E W-pi 

/l?Tll 



If-nTTl 



rn I * 'rn 



X ( $ h |"/imi(^l)d« 2 m 2 (^2) 



L 



d b 



+ all other possible full contractions |$Q h ) 



L 



^ Li 



J 



+ all other possible full contractions |$o) CO nn • 
Let us show an example how the contraction of the type ^ v (• • •) \Ev can be treated: 



(67) 



I b 1 1 



l*S> 
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f c^(*i)(- • 04 c ^) - c, 1 (t 1 )(- • O^W^^i) ^ 
1 1=! ' 1 1 

. 4(*0(- ■ -)4(**K(fi) - 4(t x )(- ■ •)c^(^)4(t i ) 

V 1 1 1 ' 1 ! / 

/ Op.{ti) 



( ^(*i)(--o4(*o c -i( t i)(---)c, l (t l ) \ 

i i i i 

v 4(*i)(--o4(**) 4(*i)(---K(*i) , 

* i ii i / 



($Sl^ 1 (ti)r 3 ^ 1 (t 1 )(---)4 1 (^)r34,(t,) ••• W • 



(68) 



Here we used c= = — c u . The summation indices (//j, v^) can be replaced by (z/«, 
Thus Eq. fl6"7j) reads 

oo 

E 



— Z 

ft 7 m! 



oo I 

-oo Z 



m=0 

x E «h E(^il^Ul^) + (^il^Ul^i)) 

X • - • 



X E K « m E (W^lmJ^m) + W^UmJ/O) 

xm h \[a hmi {ti)a l2m2 {t 2 ) ■■■ a lmmm {t m ) 
I i 1 I i 1 



+all other possible full contractions ]|$o h ) 

x(^|[*t l (i 1 )T 3 ^ 1 (tl) ••• *lJt m ) 7S*^(0^(*)*1 
I I I I I 1 =d _J | 

+ all other possible full contractions between \l/ and ^]|$o)coiin • (69) 

Now let us introduce a permutation V of the integers (1,2, • • • m,0) — > (p(l),p(2), 
, p(m), p(0)) in order to express the contractions explicitly: 



V 



( 1 2 ••• m \ 
^p(l) p(2) ••• p(m) p(0) J 

/X(l) X(2) ••• IK) \ 
^X(2) X(3) ■•• X(l) J 

'X(mx + 1) X(mx + 2) ••• X(m 2 ) \ 
X(mx + 2) X(mi + 3) ••• X{m 1 + 1) J" 

X(m P _i + 1) X(m P _i + 2) ■ ■ • X(m P ) 



x 



X 



X(m P _i + 2) X(m P _i + 3) 



X(m P _i + 1) 



(70) 



Here the decomposition of V into the product of the cyclic permutations can be 
obtained as follows: 
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i) Choose an arbitrary integer i e { 1,2, • • • m,0 } and define X(l) = i. 

ii) Pnt X(2) = p(X(l)), X(3) = p(X(2)), • • • . 

iii) If p(X(mi)) = X(l), then choose an integer j which is not yet used, and put 
X( mi + 1) = j. 

iv) Repeat the procedure. 

This is a well-known theorem in the permutation group theory. One of the cyclic 
permutations includes 0. We put 



X(m N -! + 1) = , 



(71) 



for later convenience. By applying the permutation, the summation of all possible 
full contractions between ^ and ^ can be written as 



I i — 



••• ^ m (t m )r 3 ^ m (t m )^(t)^ 
£j 1 1 J | I 



+ all other possible full contractions between \l> and \E^]|$o) 

= J2 A ^ ■ ■ ■ a p > 

■Pconn 



(72) 



where 

A l+1 = -Tr 

xr 3 

x • • 

xr 3 

x r 3 



( «CjV x(m . +1) Hx{ mi +2) tG 2 Vx(m t +l)^X(m l+ 2) 

\ l ou x(mi+1)f i x(mi+2) * oV x(mi+1) ]I x(mi+2) / ( tx{mi+1) ,t x{mi+2) ) 



jC 21 — iC 22 — — 

V V x{mi+2) Vx(m i+ 3) J (t x{mi+2) , t X{mi+3) ) 



iO 21 - iG 22 - 



X(m i+1 -1) f J 'X(m i+1 ) 



(m i+1 -l) Vx(m i+1 ) J (t X{mi+1 - 1} ,t X{mi+l) ) 
\ 



lGl VVx{n H+1 ) Vx(m l+ 1) tG 2 V X (m l+1 ) Vx(m i+ 1) 

jC 21 — iC 22 — — 

^K + .) »X(m i+ l) is x{mi+l) Llx { m i+ 1) J (t X{mi+l) , t X{mi+1) )\ 



where it is assumed that Ai + \ does not have the indices v = v§ and fi = p,Q. In the 
derivation we have used 



(73) 



i i 



^o^viHj (tii tj) iGQ 2 Ui ^. (tj, tj) 

iGo Vi^j (^»> tj) iG -y-fi. (tj, tj) 
iO 11 iG 12 - 

in 21 - iG 22 - - 

1^0 WW ^0 vim / {t . jtj) 



(74) 
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and 



M) l C d )=Tr « („,1). (75) 



The product of Go including the suffices 1/ and /x is a 2 by 2 matrix: 
^ = 



/ in 11 ir 12 — 



xr 3 



jf 21 — iG 22 — — 

\ o u x{mN _ i+2) fi x{mN _ i+3) o z/ x(mjv _ i+2) ii x(mjv _ i+3) / (t x(mjv _ i+2) ,t x( 



mjv_ 



X ■ 



xr 3 



xr 3 



* G oV x(mjv _ 1) /ix( miv ) ?G 2 i/ Jf(miv _ 1) fl X (m N ) 

iG 21 - iG 22 - - 

or 21 — ■in 22 — 



with 



V = Vq = V X (m N _ 1+ i) , (77) 
fj, = flQ = Hx(m N -i+l) • (78) 

For m in Eq. fl69|) there are ml identical terms which are obtained to each other 
by an exchange of the set of the summation indices and an integral variable (/ij, z/j, 
li, rrii, ti) <-> (fj,j, Uj, lj, rrij, tj). They are called topologically equivalent, otherwise 
the terms are called topologically distinct. One of the topologically equivalent terms 
has always 

{t 2 ) a h (t 4 ) • • • «/ m _ im m _i 

i i i i i i 

= (ajjo (a h ) iDi imuhm2 (t 1 ,t 2 ) {a h ) (a k ) iD° hm3MmA {t^U) ■ ■ ■ 

x (aj m -i)o Jo iA° m „ imm „ 1 ,« m m m (^m-i,t m ) , (79) 

where iDf z 2m2 (£i, £2) is the unperturbed phonon Green function. Thus Eq. (|69D 
becomes 

<GS,(*) ^ 22 -(t) 



fOO 



m=0,2 
X 



E (f) I dt x ---j dt v 



oc 



X E ^? 1 mi,«2m 2 ( i 1^2)^? 3m3 ,/ 4m4 (^3, i 4) * ' * « ^Lim™-l,!™mra (*m-l j O 

conn 

topologically distinct 

xA x A 2 ---A P . (80) 
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One may put the unperturbed Green functions to be diagonal: 



^Gq jy^t(t, t ) — ^Crg „(t, t , 

iGl 2 u-p(t,t') = iGl 2 u(t,t')S^ u , 

^Gq l>jl(t, t ) = ^Gg U u(t, t )5^ v , 

^Gq 77^(£, t) = iGq t )5 fMU , 



(81) 
(82) 
(83) 
(84) 
(85) 



86) 



That is, we can put 

MX(i+l) = ^X(i) , 

in Eqs. ( |73|) and (|76|). Consequently the final form of the perturbed Green function 
in the time representation is given by 

iGH(t) iG»(f) 



iGlUt) iGfJt) 



00 



x : t / 

-no \ '* / •/ -co J — 



dt r 



m=0,2 

x E E ••• E E' 

hmi hrri2 lm-\rn m -\ vwz—Vm 

x {a h ) K h {vp {1) \F} imi \vx) (a h ) K h (^(2)1^-^1^) 
x (^3)0^3(^(3)1^3^3) (^3)0 ^3(^(4)1^/3-^^4) 

X • • • 

X («/ m -i)o ^ m -i(^(m-l)|^„ _ 1 m m _i I ^m— 1 — m m _i I ^m/ 

X E Sterns (*3,*4) " " ■iA° m _ imm _ 1 (*m-l,* m ) 

conn 

topologically distinct 



xAx^a • • • 



57) 



where (fp(x(m N _i+i)) I = (^(m^)! is replaced by and |^x(mjv-i+i)) = ko) is re- 
placed by \v) in the matrix elements of F} m . Y^ Vxv% ... Vm does not contain vx{m N )- V{i) 
denotes p~ l (i). A4 takes the form 



.4 



-Tr 



r 3 



?'G 12 

^(nn-i+l) ^Xfmi-x+l) 



■^2i_ in 22 — 

V ^0 Z/ X(mi _ 1+1) Z/ X(mi l+1) ^0 Z^( m< _ 1+ 1) 



X(mj_i+1)) 



xr 3 



x 



iG 12 

^X( mi _!+2) z/ X(m l _ 1 +2) 



•^21_ iT' 22 — 
V Ztj ^X( m! _ 1+ 2) ^(m^+2) lOr ^x(m i _ 1 +2) 



X(m<_i+2)) *X(m, 



X T 3 



?'G 12 



•^21_ 1 T y 22_ 

V zCj o i/ x(m . } v x{mi) ^0 z/ x(mi) 



(il(ra,)i*IK-i+l))- 
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.4 



N 



xr 3 



V V V 

vv ^ou / (t,t x(mN _ 1+2) ) 

-p\2 _ 

<Lj o v x{mN _ i+2) i^o u x{mN _ i+2) z/ x(mjv _ 1+2) 

iG 21 - iG 22 - 

Vx imN - 1+ 2) ^X ( m N ^ 1+ 2) U x{mN _ 1+2) J (t x{mN _ 1+2h t x{mN _ 1+3) ) 



X ■ 



xr 3 



xr 3 



•/mil -/^12 _ 

<tj ^(mjv-1) <tj ^x(mjv-l) U X(m N -l) 

iG 21 - iG 22 - 

U ^ (mjv -l) u X(m N -l) U ° Vx(m N -l) ) (tx(m N -l),tx(m N )) 

or 21 — ir 22 — 
**0HH ^o/i /(t x(mjv) ,0) 



19) 



The perturbed Green function in Lehmann representation is obtained 

iG*( w ) 

lim-/ ^ e ^+«)-i)W*' ^ W "mV 

: ^+0 ft, V-oo 



5 



^H(t) iG 2 \{t) 



= E hr E 



m=0,2 

' ' conn 



/i / J-oo 27T J-oo 27T J-oo 27T 

Pc 

topologically distinct 



x E E ■•• E E' 

Zimi £377*3 Im-ltHm-1 1,2, — ,m 

x (a/Jo Mp(l)|F/ imi |l> K) ^(p(2) 14^12) 
x K)o^a(p(3)|F/ 3m3 |3) K)o^ 3 (p(4)l4-m 3 |4) ••■ 

x ("Jm-Jo Kin-Mm - l)|F/ m _ imm _Jm - 1) (a^Jo «z m _! (p(m)|F/ m _ i _ mm _ i |m) 

X J- D h mi (-^l + W p(l))j- D lm 3 (-^3 + Wp(3)) ■ ■ ■ ^ D l m -im m ^(- W m-l + <%(m-l)) 

xA\A 2 ■ ■ -Ap 

x27rh 2 8(wi + w 2 - Wp(i) - Wp(2)) 2irh 2 5{w 3 + w 4 - - Wp^) ■ ■ ■ 

x2ttTi 2 5{w m ^i +W m - Wp-(rn-l) ~ ™p{m)) • (90) 

We have used an abbreviation \i) = |z/j). As in the time-representation (j57|), 
(p(X(m N ^ 1 + 1))| = (X mN \ is replaced by and |X(tojv-i + 1)) = |0) is re- 
placed by \u) in the matrix elements of F^. Y?x2-m ^ oes no ^ contain X{mti). 
Note wo = to. A can be written 

r z q i q 

^ -I 

X ■■■T 3 -G x{m . ) (w X{ m i )) J , (91) 
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where 



For i = N, 



GS\H Gf ViVi (w) 



G?H= • (92) 



2 1 
AiV = ■^G (w )T 3 -G x{mN _ i+2) (Wx(m N ^ 1 +2)j 



1 

X ■ ■■n-^G° x{mN) {w X (m N )) , (93) 

Gxk_ 1+ i)KK-,+i)) = G° {w ) = G u {u) , (94) 
G x(mN) (w x{ m N )) = Gj(w) . (95) 

li N = 1, one may put m N _i = 0. Now the diagram rule for {1/K)G^{uj) is clear 
from Eqs. (M) - ©: 



where 



i) Draw connected topologically-distinct diagrams with the unperturbed nucleon 
Green functions, phonon Green functions and the interaction points, and as- 
sign (l/h)G°^(w), (i/h)D^ m (w), (ai) Ki(u\Fi m \fi) for emitting a phonon with 
(lm) and («;) K/(z/|F/_ m |/i) for absorbing a phonon (Im), respectively, with 
adequate suffices /ii, l\m\ etc. The definition of "connected" is that any part 
of the diagram is connected to the other part by either nucleon or phonon 
Green function. Assign intermediate energies in such a way that the interac- 
tion points conserve the energy. 

ii) Make products of the nucleon Green functions in the order indicated by the 
diagram putting r 3 at the connecting points of the nucleon Green functions. 

iii) If the product is a closed loop, take — Tr of the product. 

iv) Take summations with respect to the intermediate single-particle states, lm 
and the integral J^ 00 duJi/2Tr with respect to the intermediate energies aVs. 

ii) and iii) are the points different from the usual rule without the anomalous Green 
function. 

A simple example of diagram is shown in Fig. A-l. 

Using the diagram rule i) - iv), one can obtain the equation 

x(aj)o /^H^/j/iXaOo Ki(u\F^_ m \u 2 ) 
xi£>g>-a*), (96) 

The self-energy is defined by removing the entering and exiting nucleon Green func- 
tion from the diagram. 

We have two comments on the above derivation. The final result of the diagram 
rule is identical to that of the original formulation [fi3 |. The difference is that we 
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V t =\L 



Figure A-l. An example of diagram. 



K m oo-oi 




M^V ^ "V erf M- ny ^ 



Figure B-l. A self-energy. 



have ty. 



)ifj v ^Q. If an extension of the Nambu-Gor'kov formulation is made for a 



two-body force dependent on the spin, derivation of the diagram rule would be more 
complicated. Second, the diagram rule derived here is mathematically equivalent to 
the ways proposed in |38|. Migdal et al. do not use the 2x2 matrix form. 



Appendix B — spherically symmetric case — 

Let us consider the self-energy given by Fig. B-l under the spherical symmetry 
of the system. 

Using the diagram rule derived in Appendix A, we can write 
= f^EE r 3 ±G„(u/)T S 

Am li'm^i 

x (a x )o «a (//nv I F lm \n m v) («a)o k x (/im^ | F A f _ m l/Ary ) ^D^ m (w - a/) . 

(97) 

For the notations used, see Appendix A. Note that \i in Appendix A corresponds 
to (/i77i M ) = {nljm)^ here, m is the z-component of the multipolarity A, which is 
assumed to have only one phonon mode in this appendix for simplicity. The vertex 
matrix element can be written 

(98) 
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It is possible to make the vertex matrix elements real. 
For the unperturbed phonon Green function, we have 

iD° Xm (t) i ($ ph | T [(cL(t) + c^(t))(4_ m + Cx ^)]|$ph) 



, t > 



where the phonon creation and annihilation operators have relations 



"Am 



(") m CA 



i 

-m ■ 



(99) 
(100) 



(101) 
(102) 



|$ p h) is a phonon vacuum state. The Green function in Lehmann representation is 
given by 



iDl m {u) £ lim / dte^lHDlJt)e 

J-oo 



-5(29(t)-l)t/n 



iD» x (u>) , 



ill 



ih 



Inserting Eqs. 



and (TUB) to Eq. (fT7|), we obtain 



(103) 
(104) 



^VVr 3 ^(^)r 3 



X- 



ft 



1 



2ft A 5 A 2j„ + 1 



rft/(r) 
dr 



-Dl(u-u') 



(105) 



Appendix C — general form of particle Green function — 

The particle Green functions in the time representation are defined 

Gf(t) = -i(* \T[ Ctimft (t) ct m J|* ), (106) 
Gf(t) = -i(* o |T[4^(0cp^]|* o ), (107) 

where j^o) is an iV-particle ground state with N being an even number. The time- 
dependent creation operator of the single-particle is defined by 

r-t (t) = P iHlt / n J p-iH't/h 

H' = H- e F N . 

H is the many-body Hamiltonian. ef is the Fermi level, and N denotes the particle- 
number operator, (/xm^) indicates the time- reversed state of (fim^). Since |^ ) is 
spherically symmetric, the Green function is a scalar. Eqs. (|106| ) and ( |107| ) yield 
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i) t > 



Gf(t) 



-i^e^-W(* |c Mm J^)(^|ct |* Q ) 



iiH < 



-iEe i(E °- EiW ^^oi4^i^)(^|c^i*o) 

i 



(108) 
(109) 



G?(t) = -i{*o\{-)c^<^W Q ) 



E 



D i(Ei-E )t/h 



(*o\cpf^\*i)(*i\4frq:\*o) 



(110) 



(111) 



is the ground-state eigenenergy of H' . {^i} indicates a complete set of a many- 
body space which consists of N ± 1,N ± 3 • • -particle states. is the eigenenergy 
of \%) for 

The Green functions in Lehmann representation read 



lim 



E 



dte^' h Gf{t) e-^W-D^ 



(tto|<wJtf f ) 



u/h+ (E - Ei)/h + ir] 



E 



(112) 



G i ( w ) 



lim 



= E 

i 

Thus it is seen that we can put 
1 



+ E 









2 


u/h+(E - Ei)/h + ir] 









2 


w/h+(Ei-Ea)/1 


i — it] 



E 

a 

E 



flYJ, "T 

— cj; 



G+ 



a; + 



^10 

G+ 



W — id 



fMl 

G+ 



00 + UJq^ 



(113) 

(114) 
(115) 



The factor e luiri ( e lujr ' ) in G}-{u) ( Gjf(u>) ) comes from the definition of the Green 
functions at t = 0: 



G^(t = 0)= ilmG^t) 

^22/ 



See section 7-2 in From the time- reversal invariance it follows that 

Ei = Ej , 

f|$ )=f^ 1 |^ ) = |*o), 

{i^» i= ... = mh=... , 



(116) 
(117) 

(118) 
(119) 
(120) 
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where T is the time-reversal operator. In addition we have 



(121) 



The negative sign is due to that |^) is an odd particle-number state. By using 
Eqs. (|TLi) - ([12ID , Eq. ([TTD becomes 



G?M = E 



u/h+ (E - Ei)/h + ir] 



+ E 



u/h+ (Ei- E )/h-i7i 



E 



(*oi4 m j^) 



+ E 



<tfo|<w|tf<) 



cu/ft + (£ - Ei)/h + irj Y W A + - ^o)A - iv 



This equation implies 



G?(-a;) = -Gf(u) . 



(122) 
(123) 



For the residues we have the relations 



G+J > 



(124) 
(125) 



The anomalous Green functions are defined as 



-i(*o|T[c Mmtt (t) Cpn -]|*o) 



(126) 
(127) 



In the same manner as discussed for Gj^(uj) and G 22 (io), it turns out that 



^ (^o|c Mm J^)(^|c^|^ ) ^ (^o|c Mm J^)*(^|cp^|^o) < 



u/h+ (E - Ei)/h + ir] 



uj/h+ (Ei - E )/h - irj 



(128) 



u/h+ (E - Ei)/h + i7] 



Thus one can put 



uj/h+(Ei - E )/h - irj 



(129) 



E 



w G » = E 



G+J 



C<J — 10, 



fia 



lo — 10, 



fia 
G+ 



+ 



+ 



-00, 



fia \ 
G+J 



OJ + W 



fia 



uo + to, 



G+ 



(130) 



(131) 



with the relations 



fia \ 



Ma(~^G+; 
^a(~ W G+) 



<>G+) : 
Rfai-^Gl 



\^G 



+ 1 ' 



? 21 */ uo 



-^a ^G 



-R 



12, 



fia 

G+ 



>12*/, ,Ma 



^G 



7312*/ ua \ 



(132) 
(133) 

(134) 
(135) 
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It is worthy to note that a relation holds regardless of the time-reversal invariance 
of |tf >: 




Let us assume that the energy of |^) is equal to that of |\Ev) with a different 
particle number incidentally. Then we can introduce the basis state which has a 
mixed particle number. In this way it is possible that -R^( a; G+) an d ^ a (~ UJ G+) are 
non-zero simultaneously without causing a problem, if the system is ideally large. 



Appendix D — formula of pole and residue — 

As is mentioned in the text, we determine the real part of the pole of the per- 
turbed nucleon Green function approximately by the condition 



detG; 1 (±Reo;^ + ) 







(137) 



where G (uj) is a matrix in which the non-hermitian components of G~ (u) were 
eliminated. The advantage of this method is that the pole search is a one-dimensional 
problem. 

Using Eq. (^) in the text, we obtain 

2 



)11 , 



( d Ti 

[dZGjH^), 

Therefore it is found that 

lim 



UJ — UJ, 



fia 

G+ 



E 



d h 



(uj - uj£ a + ) 2 (uj + UJ^ + ) 



(138) 



-~-±^ a + \dujG l ^(uj) 



(139) 



and in accordance with Eq. (|137 ) we approximate 

d h 



lim 



duJG^{uj) 



(140) 



The following equation may be more convenient in the numerical calculation: 



£h'det G:\uj) 



(141) 



where S M (o;) is defined in the same way as 
For the pairing part we have 



lim 

K 2 ^ 



d h 



£K&etG-\u) 



=±Rei 



'G+ 



(142) 
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For the imaginary part of the pole it is derived that 
Irn , g „„, R]}A±^A) = Im ^M^G + ) 



Z—' I _i_r>„, , M a ' i p„, ,t*a I , 

a ; V ±Ke - u G+ ±Ke cu G+ + lv g+ 



~ Tlmu G a + . (143) 

This approximation is based on a condition that the imaginary part of the pole is 
very small. 



Appendix E — total energy — 

Derivation of the equation of the total energy is possible in the same way as that 
for the case without the pairing. (See sections 7 (Chap. 3) and 46 (Chap. 12) in [|62| ). 
We put the Hamiltonian 

Am puray.vmv 

+ -£hn x (c{ m c Xm + ~) . (144) 

Am 



The expectation value of the particle number (N) was taken in advance. For the 
particle-phonon coupling, see Appendix A. HQ\ is the phonon energy. We consider 

c^{t) = e im '% mii e- im /\ (145) 

&xm(t) = e lHt / n a Xm e- mt ' h . (146) 

Then it follows that 

^-rrc^{t) = ejtv^t) +^«Am(t) Yivm^HmWm^CvmM ■ ( 147 ) 

01 Am vm v 

Multiplying c 1 ^ (t 1 ), t' > t, from the left, using the definition of the particle Green 
function in the time representation and putting t — 0, one finds 



t=o/ 



( E ^l C lm,C^) + ( E "Am I] | K^L \ Vm^C^ C vm „ ) . (148) 



By using the Green function in Lehmann representation, it turns out that 

r°° du 1 

>m M - , -°° " Am 

E E(-^c + Xa(-^ a + ) + e F (N) + (Ht) , (149) 



-co ^r^^" + £f(a ° + ( ^ + V2)> 
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where HQ h denotes the fourth term in Eq. Q144J) . 10 in our notation always indicates 
energy. The last equation is derived for the ground state, for which the poles with 
positive infinitesimal imaginary part have negative real part. ( KeujQ^ > ) 
It is possible to rewrite the total energy by using Dyson equation as follows: 



(H) - (Ht) 
= e F (N) - i £ ± Mm / 



s F (N) 
+ 



iY 1 - lim 
& 2 



00 duj 



2vr 



-Tr 



u 
uo 



e"^ 



go 







-e 



=o 



- h a M 



e 
duo 



Y. £«(-!?+) - « £ 5 JSsJ" £ U^m\g^)^ 



h 
1 



+fr£ 22 M^G 22 Me-^ 



(150) 



We have used the relation between Gj}(u) and G 22 (u;) (Appendix C). A reasonable 
definition of the pairing energy may be 



1 t f°° doo , 
^pair = -i V r hm / — [hE 



1 



,j)^Gf(u;)e^ + hI?}(uj)-G} 2 (uj)e-^- 



1 r 



(151) 

A factor i/2 is necessary in the diagramatic derivation. If one puts hH 
hT^{ujQ^) in addition to the conditions which are used for deriving Eq. fl49D , then 
we have 

(152) 



E ■ ~ 

-^pair — 



E(2j M + l)w^ M A M . 



It is noted that Eq. ( |152|) differs from the known pairing energy in the BCS approx- 
imation by a factor 2. This is because the particle-phonon interaction introduced in 
our calculation is not of the two-body type (see Eq. (|54|), f62fl and its chap. 12). 



Appendix F — pairing gap — 

The equation to determine the pole is given by 



= h 2 Re detG" 1 (Re 

Re hT^{RBUjQ + ) 



Re ujq^ 



e° 



-ft£ 21 (Re^ + ) 



Reu£ a + + el - Re fr£ 22 (Reu;g a 



We define functions 



1 - 



^Reft£r dd ( 



uj, 



Ha \ 
G+J 



(153) 



(154) 
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v llodd/ ua \ 



ill, 



-co. 



fia 
G+ 



[CO. 



G+) 



G+, 



+ s 



11/ 



-co, 



fia 
G+ 



(155) 
(156) 



By dividing by ( Z^iReu^) ) 2 , Eq. (|153|) can be written 



Reu;^ 



£ -»+RcnSi lcvcn (R.c^ + ) 
Z M (Reu,£f + ) 



G+> 
jcur — 



Rea;^ + 



e -°+ReRSy cvcn (Bea,^) 

_____ ___ 



Z M (Re^" + ) 



(157) 



We have used £jf(^) 



^11/ 



-a;). The diagonal elements of Eq. (|157|) has the 
structure of the BCS equation for a time-reversal-invariant state. Therefore it is 
usual to define the pairing gap by 



fia \ 



ft_j, 2 (Re_£ 



158) 



( See sections 7-2 in [|3|], 10 in g4| and g7|. ) a in Eq. ( |i~58D indicates the 
quasiparticle pole. From Eq. (|157|) it is also seen that 



Reco^? 



el + RehZl ieven (Reu£ 



+ 



^(newS?) . (159) 



In the same way the perturbed single-particle energy is given by 

1 % + Betn?™(Beu%°)) . 



Z„(Reuo^) 



(160) 
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